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We investigate to what extent the magnetic properties of deposited nanostructures can be influ- 
enced by selecting as a support different surfaces of the same substrate material. Fully relativistic 
ah initio calculations were performed for Co monolayers and adatoms on Pd(lOO), Pd(lll), and 
Pd(llO) surfaces. Changing the crystallographic orientation of the surface has a moderate effect 
on the spin magnetic moment and on the number of holes in the d band, a larger effect on the 
orbital magnetic moment but sometimes a dramatic effect on the magnetocrystalline anisotropy 
energy (MAE) and on the magnetic dipole term Ta. The dependence of Ta on the magnetization 
direction a can lead to a strong apparent anisotropy of the spin magnetic moment as deduced from 
the X-ray magnetic circular dichroism (XMCD) sum rules. For systems in which the spin-orbit 
coupling is not very strong, the Ta term can be understood as arising from the differences between 
components of the spin magnetic moment associated with different magnetic quantum numbers m. 
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I. INTRODUCTION 

The magnetic properties of surface deposited nanos- 
tructures have been in the ongoing focus of many ex- 
perimental and theoretical investigations as they often 
exhibit interesting and sometimes unexpected phenom- 
ena. One of the main features in this context is that the 
local magnetic moments and their mutual interaction as 
well as the magnetocrystalline anisotropy energy (MAE) 
are in general different and often much larger in nanos- 
tructures than in corresponding bulk systems. Various 
aspects of the magnetism of many different nanostruc- 
tures were studied in the past to identify the key fac- 
tors which could then be used to tune the properties of 
such systems in a desired way. It has been known for 
some time that one such key factor is the coordination 
number, with smaller coordination numbers generally im- 
plying larger magnetic moments ji^— However, coordina- 
tion numbers alone do not fully determine magnetism of 
nanostructures. The chemical composition can play a 
significant role as well. An Fe monolayer, for instance, 
has a larger spin magnetic moment when deposited on 
Au(lll) than when deposited on Pt(lll), whereas for a 
Co monolayer it is vice versa^ The situation is even more 
diverse for the MAE where different substrates may lead 
to different properties of systems of otherwise identical 
geometries. For example, C02 and Ni2 dimers on Pt(lll) 
have out-of-plane magnetic easy axis but the same dimers 
on Au(lll) have an in-plane magnetic easy axisi^ 

Experimental research on magnetism of nanostructures 
relics heavily on the X-ray magnetic circular dichroism 
(XMCD) sum rules.—"— The strength of these sum rules 
is that they give access to spin magnetic moments /^spin 
and orbital magnetic moments //orb separately and in a 
chemically specific wayJ^ However, the XMCD spin sum 
rule does not provide /tspin alone but only its combination 



Mspin + 7Tq,, where Tq, is the magnetic dipole term (for 
the magnetization M parallel to the a axis, a=x,y, z)^ 
For bulk systems, T„ can be usually neglected but for 
surfaces and clusters the Tq, term can have significant 
influence, as it has been demonstrated experimentallySiiS 
and theoreticallyi^"— The anisotropy of the magnetic 
dipole term was predicted on general grounds^^ and some 
estimates concerning the magnitude of this anisotropy in 
non-cubic bulk systems were given based on atomic-like 
model Hamiltoniansi^ or on ab initio calculationsi^ 

Magnetic nanostructures may be prepared by combin- 
ing and arranging different magnetic elements on differ- 
ent substrates. In this respect one can also address sur- 
faces of different crystallographic orientations. Thus, it 
is important to know how the magnetic properties can be 
controlled by selecting for the substrate crystallographi- 
cally different surfaces of the same material and whether 
one can expect different effects for complete monolayers 
and for adatoms. Connected with this is the question 
about the effects on the Tq. term, because XMCD is per- 
haps the most frequently used experimental technique in 
this field and it is desirable to know how Tq, can influ- 
ence the values of magnetic moments deduced from the 
XMCD sum rules. For planning and interpreting such 
experiments, it would be very useful not only to know 
the Tq values from ab initio calculations but also to have 
a simple intuitive interpretation of the Tq term. 

In order to learn more about this, we undertook a sys- 
tematic study of Co monolayers and adatoms on Pd(lOO), 
Pd(lll), and Pd(llO) surfaces. Fully relativistic ah ini- 
tio calculations were performed to obtain /ispin; Morb, and 
Tq for different magnetization directions. The MAE was 
determined for all these systems as well. The accuracy of 
an approximative expression for the Tq term was checked 
to see whether it captures the essential physics. It is 
shown in the following that monolayers and adatoms on 



2 



different crystallographic surfaces may have indeed quite 
different magnetic properties, especially as concerns the 
MAE. Moreover, it is also demonstrated how the depen- 
dence of the Ta term on the magnetization direction leads 
to a surprisingly strong apparent anisotropy of ^spin as 
deduced from the XMCD sum rules. 



II. METHODS 
A. Investigated systems 

We investigated Co monolayers on Pd(lOO), Pd(lll) 
and Pd(llO) and also Co adatoms on the same surfaces. 
The corresponding structure diagrams are shown in Fig.[T] 
(for adatoms, obviously only one Co atom is kept). Two 
hollow adatom positions are possible for the (111) sur- 
face, differing by the position of the adatom with respect 
to the sub-surface layer; we consider the fee position in 
this work (unless specified otherwise). 

The Pd substrate has fee structure with lattice con- 
stant a=3.89 A. To determine the distances between the 
Co atoms and the substrate, we relied in most cases on 
the "constant volume approximation" : the vertical Co- 
Pd interplanar distance zco-Pd is taken as an average 
between the interlayer distance in bulk Pd and the inter- 
layer distance in a hypothetical pseudomorphically grown 
fee Co film compressed vertically in such a way that 
the atomic volume of Co is the same as in bulk CoJ^ 
In addition we took also into account relevant experi- 
mental data and results of ab-initio geometry relaxations 
when available. For example, the constant volume ap- 
proximation yields zco-Pd=l-70 A for a Co monolayer on 
Pd(lOO) while we took zco-Pd=l-65 A instead, following 
the surface X-ray diffraction experiment of Meyerhcim et 
ali^ For the other two surfaces we used the constant vol- 
ume approximation distances, namely, zco-Pd=l-96 A for 
Co on Pd (111) and zco-Pd=1.20 A for Co on Pd(llO). 
In the case of the (111) surface we can compare our 
distance with an EXAFS-derived experimental distance 
^Co-Pd=2.02 A (Ref. fish and with an ab initio equilib- 
rium distance 2Co-Pd= 1.91 A (Ref. [Ig). It follows from 
this comparison that the constant-volumc-approximation 
leads to reasonable distances. 

Systems with interplanar distances as given above will 
be called systems with "optimized geometries" . Apart 
from that, we investigate for comparison also systems 
where the Co atoms are located in ideal positions of the 
underlying Pd lattice. For this we use the designation 
"bulk-like geometry" . The interplanar distances are sum- 
marized in Tab. HI 

For adatoms we use the same zco-Pd distances as for 
monolayers. This is a simplification because the constant 
volume approximation will work worse for adatoms than 
for monolayers. For example the ab initio zco-Pd dis- 
tance for a Co adatom on Pd(lll) is 1.66 A (Ref. [13) 
in contrast to our optimized geometry value of 1.96 A. 
However, by using identical ^Co-Pd distances for mono- 



TABLE I. Vertical distances zco-Pd between the plane con- 
taining Co atoms and plane containing Pd atoms for systems 
investigated in this study. The unit is A. 



surface 


optimized geometry 


bulk-like geometry 


(100) 


1.65 


1.95 


(111) 


1.96 


2.25 


(110) 


1.20 


1.38 



layers and adatoms, the net effect due to the change in 
Co coordination can be studied. It will be shown that 
the effect of varying the distances is in fact smaller than 
the effect of monolayer-to-adatom transition. 



B. Computational scheme 

The calculations were performed within the ab initio 
spin density functional framework, relying on the local 
spin density approximation (LSDA) with the Vosko, Wilk 
and Nusair parametrization for the exchange and corre- 
lation potentiali^i The electronic structure is described, 
including all relativistic effects, by the Dirac equa- 
tion, which is solved using the spin polarized relativis- 
tic multiple-scattering or Korringa-Kohn-Rostoker (SPR- 
KKR) Green's function formalism!^ as implemented in 
the SPR-TB-KKR codei^ The potentials were treated 
within the atomic sphere approximation (ASA) and for 
the multipole expansion of the Green's function, an an- 
gular momentum cutoff ^max=3 was used. 

The electronic structure of Co monolayers on Pd sur- 
faces was calculated by means of the tight-binding or 
screened KKR techniquei^ The substrate was modeled 
by slabs of 13-14 layers (i.e. a thickness of 17-27 A, de- 
pending on the surface orientation), the vacuum was rep- 
resented by 4-5 layers of empty sites. The adatoms were 
treated as embedded impurities: first the electronic struc- 
ture of the host system (clean surface) was calculated and 
then a Dyson equation for an embedded impurity cluster 
was solved)^ The impurity cluster contains 135 sites if 
not specified otherwise; this includes a Co atom, 50-60 
Pd atoms and the rest are empty sites. 

It should be stressed that the embedded clusters define 
the region where the electronic structure and potential 
of the host is allowed to relax due to the presence of the 
adatom and not the size of the considered system. In this 
respect the Green's function approach differs from the of- 
ten used supercell approach: there is an unperturbed host 
beyond the relaxation zone in the former approach while 
in the latter approach, the supercell is terminated either 
by vacuum or by another (interfering) relaxation zone 
pertaining to an adjacent adatom. The sizes of the em- 
bedded clusters and the sizes of the supercells thus have 
a different meaning and cannot be directly compared. 

The magnetocrystalline anisotropy energy (MAE) is 

calculated by means of the torque T.^^ which describes 
the variation of the energy if the magnetization direction 
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FIG. 1. (Color online) Structure diagrams for a Co monolayer on Pd(lOO), Pd(lll) and Pd(llO). The blue and yellow circles 
represent the Co and Pd atoms, respectively. The orientation of the x and y coordinates used throughout this paper is also 
shown. 



n is infinitesimally rotated around an axis u. For uniaxial 
systems where the total energy can be approximated by 

E{e) =Eo + K2sm^{e) + Kishi^{e) , 

the difference £^(90°) — E{0°) is equal to the torque eval- 
uated for 9 = 45°!^ The torque itself was calculated by 
relying on the magnetic force theoremi^ 

Apart from the magnetocrystalline anisotropy induced 
by the spin-orbit coupling, the magnetic easy axis is 
also determined by the so-called shape anisotropy caused 
by magnetic dipole-dipole interactions. The shape 
anisotropy energy is usually evaluated classically by a 
lattice summation over the magnetostatic energy contri- 
butions of individual magnetic moments, even though 
it can be in principle obtained ab initio via a Breit 
Hamiltonian<^ In this paper, we always deal only with 
the magnetocrystalline contribution to the magnetic 
anisotropy unless stated otherwise. 

III. RESULTS 

A. Magnetic moments and magnetocrystalline 
anisotropy 

To assess the effect of selecting different crystallo- 
graphic surfaces and of going from a monolayer to an 
adatom, we calculated magnetic moments, numbers of 
holes in the Co d band and the MAE for all these sys- 
tems. The results are summarized in Tab. |IT1 For each 
system, the data are shown first for the optimized geom- 
etry and then for the bulk-like geometry (numbers in the 
brackets). The x, y, and z superscripts in the column 
header labels indicate the direction of the magnetization 
M. 

The spin magnetic moment /Xgpin and the number of 
holes in the d band nu are shown only for Afjlz, because 
they are practically independent on the magnetization 
direction: by varying it, /Xgpin can be changed by no more 
than 0.2 % and Uh by no more than 0.1 %. On the 
other hand, for /Xorb the differences can be quite large. 
The second in-plane magnetization direction Ai'jjy was 



investigated only for the (110) surface, because there is 
only very small "intraplanar anisotropy" for the (100) 
and (111) surfaces (this issue is addressed in more detail 
in Sec. IIII C\ . For bulk hep Co we get /^spin=1.61 ^s, 
/^orb=0.08 and nft,=2.48. 

Changing the surface orientation has a moderate effect 
on /Xspin and rih- The differences in /ispin when going 
from one surface to another are at most 9 %. For Uh 
these differences are at most 5 %. However, the situation 
is quite different for /iorb where the differences are 20- 
50 %. The sensitivity in /iorb finds its counterpart in 
the sensitivity of the MAE. For example, the magnetic 
easy axis for a Co monolayer is in-plane for the (100) and 
(110) surfaces but out-of-plane for the (111) surface. For 
the adatom, the easy axis is in-plane for the (110) surface 
but out-of-plane for the (100) and (111) surfaces. So in 
this respect the choice of the crystallographic surface can 
have a dramatic influence. 

Another finding emerging from Tab. |ll] is that as 
concerns /ispin, the difference between monolayers and 
adatoms is only quantitative in most cases. A surpris- 
ingly small difference in this respect is found for the (110) 
surface. As the same Co-Pd distances have been used 
for monolayers and adatoms, one observes here the net 
effect of the change in Co coordination. For //orb, the 
difference between monolayers and adatoms is obviously 
much larger than for /igpin- For the MAE this difference 
can again be essential: The magnetic easy axis for a Co 
monolayer on Pd(lOO) is in-plane while for a Co adatom 
on the same surface it is out-of-plane. Similarly, the mag- 
netic easy axis for a monolayer on Pd(llO) is parallel to 
the y-axis while for an adatom it is parallel to the x-axis. 

Changing the distance between Co atoms and the sur- 
face clearly affects the magnetic properties (cf. the values 
with and without brackets in Tab. HI)) . However, it is note- 
worthy that the effect of geometry relaxation is smaller 
than the effect of the transition from the monolayer to 
the adatom. 

We calculated also the magnetic shape anisotropy for 
the monolayers (classically, via a lattice summation, 
taking into account also moments on Pd atoms). As 
expected, this contribution favors always an in-plane 
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TABLE II. Magnetic properties of Co monolayers and adatoms on Pd(lOO), Pd(lll), and Pd(llO). The first column specifies 
whether the values are for a monolayer or for an adatom, the second column contains spin magnetic moment for the Co atom 
for A<f||z (in units of /is), the third column contains number of holes in the d band for A<f||2;. The fourth, fifth and sixth 
columns contain orbital magnetic moments for the Co atom for Ai'liz, A^Ha;, and A<f||y, respectivelly. The last three columns 
contain the MAE between indicated magnetization directions (in meV per Co atom). Numbers without brackets stand for 
systems with optimized Co-Pd distances, numbers in brackets stand for systems with a bulk- like geometry (see Sec. Ill A")) . 
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orientation of the magnetization. For Co monolayers 
on Pd(lOO) and Pd(lll), we get - E^_^ = 

—0.1 meV. For Co monolayers on Pd(llO), there is a 
small difference regarding the x and y directions: we get 

^S-dip - 4tdip - -0.07 meV and E^^^,^ - ^f^., = 
—0.09 meV. By comparing these values with the values 
shown in Tab. |TI1 we see that the shape anisotropy en- 
ergy is smaller in magnitude than the magnetocrystalline 
anisotropy energy and thus the shape anisotropy does not 
change the orientation of the magnetic easy axis as de- 
termined by the magnetocrystalline anisotropy. 



B. Induced magnetic moments 

Palladium is not magnetic as an element but it is quite 
polarizablej22i^ Spin magnetic moments induced in the 
Pd substrate by Co monolayers and adatoms are shown 
in Tab. IIIII for all three surface orientations. In the case 
of Co monolayers, the induced moments are shown for 
the first three atomic layers of Pd below the Co layer 
[denoted as Pd(l), Pd(2) and Pd(3) in Tab. Im]- Note 
that the interlayer distances are 1.95 A, 2.25 A and 1.38 A 
for the (100), (111) and (110) surfaces, respectively. The 
relatively large /ispin for the Pd(2) and Pd(3) sites in 
the case of the (110) surface reflects the relatively small 
interlayer distance for this crystallographic orientation. 

In the case of adatoms, the description is formally more 
complicated because Pd atoms belonging to the same co- 
ordination shell around the Co atom are not all equiva- 
lent: some of them belong to the surface layer, some to 
the sub-surface layer and so on. In order not be over- 



TABLE III. Spin magnetic moments for Pd atoms which are 
first, second and third nearest neighbors of Co atoms, in units 
of fiB ■ As in Tab. |TT1 the numbers without brackets stand for 
systems with optimized geometry and the numbers in brackets 
stand for systems with bulk-like geometry. 





Pd(l) 
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Co on 


Pd(lOO) 




monolayer 


0.29 


0.17 
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monolayer 


0.29 


0.22 


0.17 




(0.29) 


(0.24) 


(0.19) 


adatom 


0.15 


0.04 


0.04 




(0.15) 


(0.05) 


(0.04) 



whelmed by too much data, we display here only mo- 
ments averaged over all atoms of a given coordination 
shell. Symbols Pd(l), Pd(2), and Pd(3) in Tab.HJstand 
now for the first, second, and third shell of Pd atoms 
around the Co adatom. 

Moreover, we also calculated the orbital magnetic mo- 
ments for the Pd atoms in all systems and we found that 
fiorh amounts to about 8-17 % of the corresponding /ispin- 
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In this section we deal only with magnetic moments 
on those Pd atoms which are close to the Co atoms. The 
issue of more distant Pd atoms and of the total charge 
contained in the polarization cloud is dealt with in the 
Appendix. Here, we would only like to stress that it fol- 
lows from the analysis outlined in the Appendix that our 
model system is clearly adequate to yield reliable values 
of induced magnetic moments for the Pd(l), Pd(2), and 
Pd(3) sites. 

C. Azimuthal dependence of MAE 

In general, the MAE defined as the difference between 
total energies for in-plane and out-of-plane orientation of 
the magnetization will depend on the azimuthal angle (j). 
This dependence is often ignored but may sometimes be 
significant. In our case, the intraplanar MAE E'^^^ — i?^^^ 
is quite comparable to E'^^^ — E^'^^ or E^^^ — E^^^ for 
the (110) surface (see Tab. |lT]). To get a more com- 
plete picture, we inspect the azimuthal dependence of 
i;(ll)((/)) - where E^^^^c/)) is the total energy if M 

is in the surface plane (9=0°) with the azimuthal angle 
(j>. Our results for a Co adatom on all three Pd surfaces 
are shown in Fig. [51 The data reported here were ob- 
tained for the bulk-like geometry but the trends would 
be similar for any zco-Pd distance. 

One can see from Fig.[5]that the E^^^'>{(f>) — E'^^^ curves 
follow the symmetry of the appropriate surface, as ex- 
pected. The amplitude of these curves is the most in- 
teresting information here. For high-symmetry surfaces, 
it is almost negligible: 0.008 meV or 3 % of the average 
value for Co on Pd(lOO) and 0.06 meV or 1 % of the 
average value for Co on Pd(lll). For the (110) surface, 
however, the amplitude is 2.6 meV and to speak about an 
average MAE does not make sense in this illus- 
trated by the fact that the magnetic easy axis is in-plane 
for (/) = 0° and out-of-plane for (f) = 90°. 

D. Relation between magnetic dipole term and 
m-decomposed spin magnetic moment 

The spin magnetic moment sum rule for the ^2,3 edge 
XMCD spectra can be written for a sample magnetized 
along the a direction as^ 

I / (A/.i3 - 2A^.,J AE = + , (1) 
-f J nil 

where A/j.^^ , are the differences A/x = — be- 
tween absorption coefficients for the left and right circu- 



larly polarized light propagating along the a direction, / 
is the integrated isotropic absorption spectrum, /^spin is 
the local spin magnetic moment (only its d component 
enters here), Uh is the number of holes in the d band, and 
Tq, is the magnetic dipole term related to the d electrons. 
Ta can be written as^i^ 

T - if ) 

— ^ \-'-a/ 1 




where 

Qc.,3 - - Sr^rO (3) 

is the quadrupole moment operator and Sa is the spin 
operator. If z is the quantization axis, the eigenvalues of 
S;, are ±(l/2);i. 

A more transparent expression for Ta can be obtained 
if the spin-orbit coupling can be neglected. Then one can 
write^ 

t. = (-^) Qx.S, for M\\x , 

fy = {-^)QyyS, for M\\y , (4) 

t = [-^)q.J-. for M\\z , 

where Qxx, Qyy and Qzz are quadrupole moment com- 
ponents referred to the crystal (global) reference frame 
and Sz is the spin component with respect to the local 
reference frame in which z is identical to the spin quan- 
tization axis. We are interested in the expectation value 
of the Ta operator acting on the d components of the 
wave function in the vicinity of the photoabsorbing site. 
Using for the sake of clarity a simplified two-component 
formulation instead of the full Dirac approach, the wave 
function can be expanded in the angular-momentum ba- 
sis as 

^Ekir) = J2J2^ELnir)YUr)x^'^ (5) 

£m s 

to obtain 

Ta = (-^) f 'dE f dk {i:Ek\QaJz\^Ek) ■ (6) 

Restricting ourselves just to the £ = 2 component and 
omitting the corresponding subscript in a^^^^(r), we get 
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FIG. 2. (Color online) Difference between total energies for in-plane and out-of-plane magnetization for a Co adatom on Pd 
(100), (111), and (110) surfaces (bulk-like geometry). Points are results of the calculation, dashed lines are sinusoidal fits. The 
orientation of the x and y axes is as in Fig. [1] 
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Ep 
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(^2?7i I Qaa I ^2m') 



(7) 



where the spin-dependent number of states matrix N^^^, 
is defined as 



(^) 



BZ 



dE dk rMra'^llir)alXEir) 



The difference of the diagonal terms of NmL is just 
the spin magnetic moment decomposed according to the 
magnetic quantum number to, 

(^A'-B) {-^mm ^ ^mm) ~ Mspin ' 

with the sum of all the m components giving the total 
spin magnetic moment (of the d electrons, in our case). 
Therefore, if it was possible to restrict the sum ([7]) just 
to the terms diagonal in m, one would have 



(8) 



The procedure we employed above is essentially the one 
suggested by Stohr,— but we present it here in a more 
explicit way. 



TABLE IV. Diagonal components of the quadrupole opera- 
tor in the basis of real spherical harmonics. Non-diagonal 
components are all zero except for the components given in 
Eq. (O. 



0^yz\Qaa l^s) 
{V322_r2 \Qaa\Y^^-2_^2) 
(Yxz \ Qaa I "^^2 ) 



Qxx 
_ 2 
7 

4 
7 
2 
7 

_2 

7 
2 



Qyy 
_ 2 
7 

_ 2 
7 
2 
7 
4 
7 

_ 2 
7 



The coefficients {Y2m\Qaa\Y2m') can be obtained by 
analytic integration. If we use the basis of real spherical 
harmonics, the only "cross-terms" which are non-zero are 



{Yx2-y2\Qxx\Y3^2_, 



{Y,:2^y2\Qyy\Y,,2_,2) = -i2/7)V3 



(9) 



Otherwise, only the diagonal terms {Y2m\Qaa\Y2m) are 
non-zero and we list them in Tab. IIVI (see also Refs. [l3 
andlsil). Therefore, in the absence of spin- orbit coupling, 
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Eq. ([5]) presents an exact expression for and an ap- 
proximate expression for and Ty [due to the existence 
of non-diagonal terms Q]. As argued by Stohr;^ the 
non-diagonal terms drop out of the sum in Eq. ([7|) for 
high symmetry systems. 

Eq. ([HI) together with Tab. IIVI illustrate the common 
statement that the magnetic dipole term is related 
to spin anisotropy: if the m-components of /igpin are all 
identical, is zero (in the absence of spin-orbit cou- 
pling). It is also evident from Eq. ([8|) and Tab. IIVI that 
the Ta term will generally depend on the magnetization 
direction a. 

To get a more quantitative feeling of how the various 
contributions add together to generate Ta , we present in 
Tab. |V] the m-decomposed magnetic moment /x^™^ and 
individual terms of the sum ([8|) for Co monolayers on 
Pd surfaces. One can see that the term is formed by 
a competition between those m components which con- 
tain the a coordinate and those which do not (they con- 
tribute with an opposite sign, as it can be seen also from 
Tab. II Vp . In fact, this is what is meant by the statement 
that the Ta term describes the anisotropy of ^spin- 

Eq. (O gives an intuitive insight into Ta provided that 
the underlying approximations — the neglect of the spin- 
orbit coupling and of the non-diagonal terms shown in 
Eq. ([9]) — are not too crude. To check this, we com- 
pare the values of Ta calculated via the exact relation in 
Eq. ^ and via the approximative Eq. (jS]). Special at- 
tention is paid to the differences between the Ta terms 
for different orientations of M, because the 7{Ta — Tp) 
quantities determine the apparent anisotropy of /igpin as 
deduced from the XMCD sum rule in Eq. ([!]). The out- 
come for both monolayers and adatoms is summarized 
in Tab. IVII Let us recall that for bulk hep Co, the mag- 
netic dipole term is very small (we get T^^— 0.002 /is). 
Note that all values presented in Tabs. IVtiVIl were ob- 
tained from fully relativistic calculations, including the 
spin-orbit coupling. 

One can see from our results that the approximative 
expression for Ta works quite well for the Co-Pd systems: 
quantitative deviations sometimes occur but the main 
trend is well maintained. One can expect that for systems 
with a strong spin-orbit coupling the deviations between 
Eqs. (0) and ® will be larger. 

The last two columns of Tab. IVII contain the values of 
7(Ta; — T^) and, for the case of the (110) surface, also of 
7(Ty — Tz). These values are comparable to /Xspin which 
means that even though /igpin practically does not depend 
on the magnetization direction at all, its combination 
Mspin + 7Ta probed by the XMCD sum rule may strongly 
depend on the magnetization direction. 



IV. DISCUSSION 

We investigated how the magnetic properties of Co 
adatoms and monolayers can be manipulated by select- 
ing different supporting Pd surfaces. We found that this 



has a moderate effect on ^spin and Uh, larger effect on 
/iorb and dramatic effect on the MAE and on the Ta 
term. For the adatoms the effect is larger than for the 
monolayers. Moreover, the transition from monolayers 
to adatoms has a larger effect than a moderate variation 
in the height of the Co layer above the substrate. If the 
spin-orbit coupling is not very strong, the Ta term can be 
understood as arising from a competition between those 
TO-decomposed components of ^spin which are associated 
with the a coordinate and those which arc not. 

In the past, the influence of the orientation of super- 
lattices (multilayers) on magnetic properties was already 
investigated, however, the focus was mainly on the role 
of defects and interface abruptnessi^ Here, we deal with 
perfect monolayers and surfaces and investigate how sole 
selection of a different surface can affect various quanti- 
ties related to magnetism. Likewise, the importance of 
the Tz term for an XMCD sum rules analysis has been 
highlighted before when it was found that the absolute 
value of 7Tz amounts to about 20 % of /Xspin for some 
low-dimensional systemsi^ or that for atomic clusters 
Mspin can show a different behavior with changing clus- 
ter size when compared to //spin + 7Tz In this study 
the importance of the anisotropy of the magnetic dipole 
term in nanostructures is stressed for the first time and 
it should be noted that the anisotropy of Ta which we 
highlight here is primarily connected with the breaking 
of the crystal symmetry at the surface and occurs even 
without spin-orbit coupling. 

For the monolayers, the changes in //spin when going 
from one surface to another reflect the corresponding 
changes in the coordination numbers: /tspin is largest for 
the (110) monolayer where each Co atom has got only 
two nearest neighboring Co atoms, next comes the (100) 
monolayer with four Co neighbors and the lowest fj-spin 
is obtained for the (111) monolayer with six Co neigh- 
bors. This complements an analogous trend found ear- 
lier for freoi and supported clustersj^i^i^ The magnetic 
moments induced at individual Pd atoms are larger for 
Co monolayers than for Co adatoms, which reflects the 
fact that for monolayers, Pd atoms are polarized by more 
than one Co atom. 

The large amount of data gathered here for quite a 
complete set of systems allows a comprehensive look at 
the relation between the MAE and the anisotropy of /(orb- 
In this respect Bruno's formula^^ 



[a) 
Morb 



/^orb 



(10) 



connecting the differences of total energies to the dif- 
ferences of orbital magnetic moments for two orienta- 
tions of the magnetization, a and /3, proved to be very 
useful2& despite its limitations)^ which become more se- 
vere in the case of multicomponent systems with large 
spin-orbit coupling parameter ^ for the non-magnetic 
componenti^i^ To assess the situation for 3d-Ad alloys, 
we compare the differences A/Xorb and A.E, using all 
the appropriate values given in Tab. |lTl The outcome is 
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TABLE V. Spin magnetic moment decomposed according to the magnetic quantum number m together with the corresponding 
Tq""' = ^nl^i^{Y2m\Qaa\y2m) terms of the decomposition ((S)) for Co monolayers on Pd (optimized geometry). The sums of 
these components are shown in the last row for each system and they correspond to the total ^spin, Tz, T^, and Ty of the d 
electrons [evaluated using the approximative expression (|8]l in the case of Tq]. 

component fij^^j rj'"' tJ'"' T^'"' 

Co on Pd(lOO) 

xy 0.319 0.092 -0.046 -0.046 

yz 0.465 -0.066 0.133 -0.066 

3z^-r'^ 0.365 -0.104 0.052 0.052 

xz 0.465 -0.066 -0.066 0.133 

x'^-y^ 0.449 0.128 -0.064 -0.064 

sum 2.062 -0.018 0.009 0.009 

Co on Pd(lll) 

xy 0.339 0.097 -0.048 -0.048 

yz 0.428 -0.061 0.122 -0.061 

3z^-r'^ 0.490 -0.140 0.070 0.070 

xz 0.428 -0.061 -0.061 0.122 

x^ - y^ 0.339 0.097 -0.048 -0.048 

sum 2.023 -0.069 0.034 0.034 

Co on Pd(llO) 

xy 0.397 0.113 -0.057 -0.057 

yz 0.346 -0.049 0.099 -0.049 

3z^-r'^ 0.515 -0.147 0.074 0.074 

xz 0.527 -0.075 -0.075 0.151 

x'^ - 0.343 0.098 -0.049 -0.049 

sum 2.128 -0.060 -0.009 0.069 



TABLE VL Magnetic dipole term for Co monolayers and adatoms on Pd(lOO), Pd(lll) and Pd(llO) (optimized geometries) 
for different magnetization directions. For each system, the first line ("exact") contains values calculated using Eq. ((2| and the 
second line ("approx.") contains values calculated using Eq. ([S|. The Ty terms were evaluated only for the (110) surface. 











Ty 


7(r, -T,) 


7{Ty-n 








Co on Pd(lOO) 








monolayer 


exact 


-0.017 


0.010 




0.188 






approx. 


-0.018 


0.009 




0.184 




adatom 


exact 


-0.024 


0.015 




0.275 






approx. 


-0.026 


0.013 
Co on Pd(lll) 




0.276 




monolayer 


exact 


-0.066 


0.035 




0.707 






approx. 


-0.069 


0.034 




0.723 




adatom 


exact 


-0.146 


0.080 




1.577 






approx. 


-0.154 


0.077 

Co on Pd(llO) 




1.618 




monolayer 


exact 


-0.057 


-0.008 


0.068 


0.339 


0.872 




approx. 


-0.060 


-0.009 


0.069 


0.360 


0.904 


adatom 


exact 


-0.112 


-0.020 


0.141 


0.644 


1.768 




approx. 


-0.117 


0.011 


0.106 


0.900 


1.566 



shown in Fig. [31 together with a straight line represent- 
ing Eq. (jlO[) . Here we take 85 meV for the spin-orbit 
coupling parameter ^ [which appears to be a rather uni- 
versal value for Co as our calculations yield ^ of 85.4 mcV, 
84.5 meV, 84.9 meV and 85.1 meV for bulk hep Co and 
for a Co monolayer on Pd(lOO), Pd(lll) and Pd(llO), 
respectively]. It follows from Fig. [3] that Bruno's for- 
mula Eq. pH)) works quite well for adatoms (albeit with 



some "noise") but not so well for monolayers, where re- 
lying solely on Eq. (jlOp might even lead to a wrong sign 
of the MAE. This may be connected with the fact that 
for monolayers, the MAE is generally not very large and 
hence small absolute deviations from the rule given in 
Eq. (fTO)) can lead to large relative errors. 

The sizable intraplanar anisotropy i?'-^-' — which 
we get for a Co monolayer on Pd(llO) had to be expected 
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FIG. 3. (Color online) Dependence of the MAE for Co mono- 
layers and adatoms on the difference of orbital magnetic mo- 
ments for respective magnetization directions. The dashed 
line represents Bruno's formula in Eq. (jlOp . 



as this system could be viewed as a set of Co wires which 
are surely anisotropic in this respect. However, we get a 
very strong azimuthal dependence of the MAE also for 
the adatom on the (110) surface which is quite surprising 
as this can be only caused by the underlying substrate. 
The magnetic moments at Pd atoms are not very large 
(Tab. IIII[) . neither is the spin-orbit coupling parameter 
^ for Pd in comparison to, say, 5d elements. Thus, this 
seems to be yet another example of the extreme sensi- 
tivity of the MAE. At the same time, let us note that 
the calculated azimuthal dependence of the MAE can be 
accurately fitted by smooth sinusoidal curves (see Fig. 
(2) which indicates a very good numerical stability of the 
computational procedure. 

The intraplanar anisotropy for a Co adatom on the 
Pd(lll) surface can be compared to similar systems in- 
vestigated in the past. In particular, for a Co adatom 
on Pt(lll) the amplitude of the i?(ll)(0) — E^^^ curve is 
about 2 % of the average value^^ i.e., similar to the cur- 
rent case. For a 2x2 surface supercell coverage of Fe on 
Pt(lll), this amplitude is 10-25 % (depending on the 
geometry relaxation)^ but this situation is already quite 
distinct from the isolated adatom case. 

According to our calculations, a Co monolayer on 
Pd(lOO) has an in-plane magnetic easy axis, a Co mono- 
layer on Pd(lll) has an out-of-plane magnetic easy axis 
and the difference between the respective MAE values is 
about 1 meV, which can be seen as a measure of how 
much the out-of-plane magnetization is preferred by the 
Co/Pd(lll) system in comparison with the Co/Pd(100) 
system. This is similar to what was calculated for Co/Pd 
multilayers: both CoiPdg (100) and CoiPda (111) mul- 
tilayers have an out-of-plane magnetic easy axis but the 



MAE per unit cell is by about 0.9 meV larger for the 
(111) multilayer than for the (100) multilayer 4^ 

The theoretical values for the anisotropy of Tq. shown 
in Tab. IVII can be compared with experimental data for 
a similar system, namely, a single Co(lll) layer sand- 
wiched between two thick Au layers. By extrapolat- 
ing results obtained via angle-dependent XMCD mea- 
surements, Weller et alM- obtained 7Tx = 0.43 fiB and 
7Tz = —0.86 i-iB- Our values for a Co monolayer on 
Pd(lll), 7T^ = 0.24 i^iB and 7T, = -0.46 fis (see Tab. 
IVip . are fully consistent with this. 

We expect that our values for /iorb will be systemat- 
ically smaller than experimental values because we rely 
in the LSD A which usually underestimates A'orb'^'^ The 
same may be also true for the MAE. However, this does 
not affect our conclusions. 

We used potentials subject to the ASA which may limit 
the numerical accuracy of our results, particularly as con- 
cerns the MAE. On the other hand, our results do not 
differ too much from results of full-potential calculations, 
especially in the case of monolayers. For a Co monolayer 
on Pd(lOO), we get an in-plane magnetic easy axis with 
an MAE of -0.73 meV per Co atom while Wu et al^ ob- 
tained for the same zco-Pd distance (1.65 A) a theoretical 
MAE of -0.75 meV. Magneto-optic Kerr measurementsi^ 
as well as XMCD experiments^ showed that the mag- 
netic easy axis of ultrathin Co films on Pd(lOO) is indeed 
in-plane (the experiment includes also an in-plane con- 
tribution from the shape anisotropy) . Note that the the- 
oretical MAE of -0.18 meV given in Ref. [13 was obtained 
for a partially disordered Co monolayer simulating the 
growth conditions, so it cannot be directly compared to 
our results obtained for an ideal monolayer. 

For a Co monolayer on Pd(lll), we get a ^spin value of 
2.01 fiB ill a Co ASA sphere with a radius of 1.46 A while 
the full-potential calculations of Wu et alM- led to a /igpin 
value of 1.88 /is obtained within a Co muffin-tin sphere 
with a radius of 1.06 A — both calculations thus again 
give consistent results. For Pd atoms just below the Co 
layer, we get a /Xspm value of 0.32 /is in a sphere with a 
radius of 1.49 A while the corresponding /igpin value of 
Wu et al^ obtained within a sphere having a radius of 
1.32 A is 0.37 ^B- In this last case, one has to bear in 
mind that Wu et al^ used a thin slab of only five Pd 
layers sandwiched between two Co layers which clearly 
favors a larger Pd polarization in comparison with just a 
single Co-Pd interface considered in this work. 

For adatoms, the ASA may be more severe than for 
monolayers, nevertheless, the agreement between our cal- 
culations and the results obtained via a full potential cal- 
culation is pretty good (see the end of the Appendix) . As 
a whole, the accuracy of our calculations is sufficient to 
warrant the conclusions which rely on comparing a large 
set of data and not only on results for a singular system. 

It follows from our results that one can change the 
magnetic easy axis from in-plane to out-of-plane direction 
just by using as a substrate another surface of the same 
element. This could be used as yet another ingredient 
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for engineering the MAE of nanostructures, which has 
become a great challenge recently 4^ We also showed that 
the magnetic dipole Tq term can mimic a large anisotropy 
of /^spin as determined from the XMCD sum rules. Hence, 
the anisotropy of has to be taken fully into account 
when analyzing XMCD experiments on nanostructures. 

V. CONCLUSIONS 

Co monolayers and adatoms adsorbed on different sur- 
faces of Pd exhibit quite different magnetic properties. 
The effect on /igpin is moderate, the effect on /ioib is larger 
while the effect on the MAE and on the magnetic dipole 
term may be crucial. A surprisingly strong azimuthal 
dependence of the MAE is predicted for a Co adatom on 
Pd(llO). 

The dependence of on the direction of the magne- 
tization can lead to an apparent anisotropy of the spin 
magnetic moment as deduced from the XMCD sum rules. 
For systems with small spin-orbit coupling, the Ta term 
can be related to the differences between components of 
the spin magnetic moment associated with different mag- 
netic quantum numbers. 



T ' I ' I ' I ' I ' I ' I ^ 




2 4 6 8 10 12 

radius of zone [A] 

FIG. 4. (Color online) Sum of the spin magnetic moments at 
the Co adatom and at those substrate Pd atoms which are 
enclosed in hemispherical zones of the given radii, for four 
embedded cluster sizes (identified by numbers of Pd atoms 
contained in them). 
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Appendix: Effect of the size of the relaxation zone 

When studying the magnetism of adatoms, one should 
address the question to which extent the host around the 
adatom has to be allowed to polarize. Zeller showed^ 
that the polarization cloud around a magnetic impurity 
in bulk Pd extends at least up to 1000 atoms. Sipr et al^ 
showed that the convergence of the MAE with respect to 
the slab thickness and/or with respect to the size of the 
supercell which simulates the adatom is much slower than 
the convergence of magnetic moments. In view of these 
facts, it is desirable to explore more deeply the situation 
for the systems considered in this work. 

As a test case, we select a Co adatom on Pd(lll). 
To facilitate the comparison with calculations done by 
other methods, we put the Co adatom in an hep hollow 
site, with the vertical distance between the Co adatom 
and the Pd surface layer as zco-Pd=l-64 A. Our system 
is thus similar to the system investigated by Blohski et 
al^ (the main difference with respect to Ref. [13 is that 
we do not consider any buckling of the substrate). To 



check the convergence with respect to the size of the zone 
where the electronic structure is relaxed, we probed a 
series of embedded cluster sizes, starting with relaxing 
the electronic structure just in three Pd atoms (i.e., up 
to the distance of 2.3 A from the Co adatom) and ending 
with relaxing it in 220 Pd atoms (up to 11.7 A from the 
Co adatom) . To safely accommodate this large embedded 
clusters, we model the Pd substrate by a slab of 19 layers 
[contrary to 13 layers used in other calculations involving 
the Pd(lll) surface in this work]. The largest embedded 
cluster with 220 Pd atoms contains Pd atoms located 
within the fith layer below the surface and comprises 329 
sites altogether. 

First we investigate the convergence of the spin mag- 
netic moments. This can be achieved by inspecting the 
total /ispin contained inside a hemisphere stretching from 
the adatom up to a certain radius. The dependence of 
this total /igpin on the radius of the hemisphere forms an 
"integral magnetic profile" . This is presented in Fig. 2] 
for four embedded cluster sizes containing 7, 46, 133, and 
220 Pd atoms, respectively. The total /.tspin for a sphere 
with zero radius is obviously just the /igpin value of the 
Co adatom. With increasing sphere radius the spin mag- 
netic moments of enclosed Pd atoms are added to it. If 
the radius of the hemisphere becomes larger than the ra- 
dius of the embedded cluster, the total /Xgpin obviously 
does not change any more because the Pd atoms outside 
the embedded impurity cluster are nonmagnetic. 

It follows from Fig. 2] that the spin magnetic moment 
of the adatom as well as magnetic moments induced in 
the nearest Pd atoms are actually already well described 
by relatively small embedded clusters. However, the total 
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50 100 150 200 

Pd atoms in embedded cluster 

FIG. 5. The MAE of a Co adatom in an hep position on 
Pd(lll) for different sizes of the embedded clusters. 



Mspin converges only very slowly with increasing size of 
the relaxation zone because even quite distant Pd atoms 
still contribute with their non-zero //spin- Our results sug- 
gest that the magnetic moments on all the Pd atoms do 
not arise due to a direct interaction with the Co adatom. 
Rather, the adatom induces a magnetization in its near- 
est neighbors, then these further induce magnetization in 
the next coordination shell and so on. The emerging pic- 
ture of how the magnetism spreads through the Pd host 
is thus consistent with the picture suggested by Polesya 
et ali^ in terms of an exchange-enhanced magnetic sus- 
ceptibility (see Fig. 4 of Ref. Isol and the associated text). 
A plot analogous to Fig. |4] could also be drawn for /iorb 
exhibiting the same features as seen in Fig. 21 

Our results on the convergence of the magnetic mo- 
ments may raise objections about the convergence of the 
MAE. If embedded clusters containing as much as 220 



Pd atoms still do not fully account for the host polariza- 
tion, can one get reliable results for the MAE, which is 
sensitive to the way the substrate is treated?— To check 
this, we calculated the MAE for a series of embedded 
cluster sizes (Fig. [5]) . One can see that in fact the MAE 
converges quickly with increasing size of the embedded 
cluster. Already with a relaxation zone including only 
46 Pd atoms, which corresponds to a radius of the hemi- 
sphere of 6.9 A containing Pd from up to the third Pd 
layer below the surface, the accuracy of the MAE is bet- 
ter than 1 %. This means that all the results presented 
in this work are well converged. 

The data in Fig. [5] demonstrate that it is sufficient 
to include a rather small polarization cloud within the 
Pd host in order to get convergence in the MAE values. 
More distant Pd atoms do not contribute to the MAE, 
even if they are magnetically polarized. This conclusion 
is not in contradiction with an earlier result that reliable 
values of the MAE can be obtained only if the host is 
represented by slabs of at least ten layers^ because that 
result concerned the total "physical" size of the model 
system while in this appendix we focus only on the size 
of the zone where the electronic structure is allowed to 
relax to the presence of an adatom (or of an adsorbed 
monolayer). 

To complete this part, we should compare our results 
with the results of Blonski et al^ which were obtained 
by performing a plane-wave projector-augmented wave 
(PAW) calculation for a supercell comprising five-layers 
thick slabs and a 5x5 surface unit cell. As concerns the 
Co adatom itself, /ispin and /iorb for the in-plane magneti- 
zation direction and /Xorb for the out-of-plane magnetiza- 
tion direction are 2.48 /is, 0.15 /is, and 0.27 /is in this 
work and 2.24 fj,B, 0.19 /ts, and 0.22 in Blonski et 
al^ As concerns the MAE calculated via the magnetic 
force theorem (torque method), it is 1.90 meV out-of- 
plane in this work and 0.72 meV out-of-plane in Blonski 
et alM' The value for the induced /tspin in the nearest Pd 
atoms is 0.28 jis in this work and 0.33 /is in Bloiiski et 
al^ All these values are in good agreement, considering 
the differences between both approaches. 
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